Spin cloud induced around an elastic scatterer by the intrinsic spin-Hall effect 
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Similar to the Landauer electric dipole created around an impurity by the electric current, a spin 
polarized cloud of electrons can be induced by the intrinsic spin-Hall effect near a spin independent 
elastic scatterer. It is shown that in the ballistic range around the impurity, such a cloud appears 
in the case of Rashba spin-orbit interaction, even though the bulk spin-Hall current is absent. 
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The spin-Hall effect attracts much interest because it 
provides a method for manipulating electron spins by 
electric gates, incorporating thus spin transport into con- 
ventional semiconductor electronics. As it has been ini- 
tially predicted, the electric field E induces the spin flux 
of electrons or holes flowing in the direction perpendicu- 
lar to E. This spin flux can be due either to the intrinsic 
spin-orbit interaction (SOI) inherent to a crystalline solid 
|lj|, or to spin dependent scattering from impurities [2|- 
Intrinsic spin-Hall effect corresponding to the former sit- 
uation has been observed in p-doped 2D semiconductor 
quantum wells jjj] , while the extrinsic effect related to the 
latter scenario has been detected in n-doped 3D semicon- 
ductor films [2|. 

Most of the theoretical studies on the spin-Hall effect 
(SHE) has been focused on calculation of the spin current 
(for a review see p|). On the other hand, since the spin 
current carries the spin polarization, one would expect a 
buildup of the spin density near the sample boundaries. 
Such a spin accumulation near interfaces of various na- 
ture was calculated in a number of works |y, Li l3 • This 
accumulated polarization is a first evidence of SHE that 
has been observed experimentally in Ref. [ilJI. In fact, 
measuring spin polarization is thus far the only practical 
way to detect SHE. 

Yet the spin accumulation near interfaces is not the 
only signature of SHE. To draw an analogy with the 
charge transport, one can expect that similar to Lan- 
dauer charge dipoles created by the DC electric current 
around impurities |2|, nonequilibrium spin dipoles must 
be formed subsequent to the spin-Hall current. One may 
expect that the spin cloud will appear around a spin-orbit 
scatterer in case of extrinsic SHE, as well as around a 
spin-independent scatterer, in case of the intrinsic effect. 
We will consider the latter possibility for a 2D electron 
gas with Rashba interaction. The polarization in the 
direction perpendicular to 2DEG will be calculated in 
the ballistic range around an impurity represented by an 
isotropic spin independent scattering potential. Besides 
conventional semiconductor quantum wells this analysis 
can be applied to metal adsorbate systems with strong 
Rashba type spin splitting of surface states [l3| . In this 
case the spin cloud can be measured by STM with a mag- 



netic tip. 

The Landauer electric dipole has been calculated 
jld lllj basing on the asymptotic expansion of the elec- 
tron waves elastically scattered from an isolated impurity. 
Subsequent averaging of the corresponding spatial prob- 
ability weighted by the Boltzmann distribution function 
of incident wavevectors produces the dipole distribution. 
The spin cloud could be obtained in a similar way. In- 
stead, we choose a Green function method combined with 
the linear response theory. Within this method the spin 
density is given by the standard Kubo formula where the 
scattering potential of a target impurity, at a fixed po- 
sition Vi, is included into the Green functions, up to the 
second perturbation order. Other impurities are assumed 
to be randomly distributed over a 2D plane, so that the 
calculated spin density is averaged over their positions. 

We assume that a uniform external electric field is ap- 
plied to 2DEG. The field is represented by the vector 
potential A, E = iwA/c, with oj ^ in the DC regime. 
The corresponding interaction Hamiltonian is eA • v/c, 
where the velocity v^ , j — x,y, includes the spin-orbit 
correction 9(hk • cr^jdW . The spin-orbit field hk is a 
function of the two-dimensional wave-vector k. In its 
turn, the spin-orbit interaction is written in the form 



B, 



hk • cr, 



(1) 



where (T=(cr^, a-y ,ct^) is the Pauli matrix vector. We as- 
sume that the target impurity, located at r^, has a scat- 
tering potential U{r — r^). In 2D geometry the corre- 
sponding Born amplitude is given by 



/(k,k') 



dr^U{r)e 



;(k-k': 



(2) 
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where ft = 1 and ip is the angle between k and k'. Both 
the scattered and the incident wavevectors are taken at 
the Fermi circle with the radius kp- Other impurities, 
which not necessarily are of the same nature as the tar- 
get impurity, are randomly distributed within a sample. 
They create the random potential T4c(r) which is as- 
sumed to be delta correlated, so that the pair correlator 
{Vsc{r)Vscir')) = TSir - r')/7r7Vi., where Np is the elec- 
tron density of states at the Fermi energy, and F — 1/2t 
is expressed via the mean elastic scattering time r. As- 
suming that the semiclassical approximation EpT ^ 1 



is valid, one can apply the standard perturbation theory 
[la. Il4| when calculating the configurational averages of 
Green functions and their products. In the leading or- 
der of {EpT)~^ and up to the second order in t/(r — r^), 
the average retarded Green function in the momentum 
representation is given by 

GLu' (^) = 6^^^'GT (^) + GLlV (c.) + G"-^^} (u;) , (3) 
with the unperturbed function given by the 2x2 matrix 



where Ei^—k"^ /{2m*). Other functions in (Q are 



(4) 



(5) 



GLl'^ iu^) = Gl^'^ (c.) V U^^„Gl^^^ (c.)t/k"k'G;,^"^ (c.) . 



The matrix elements t/k'k = — \/27rfcF/(k, k') exp[i(k — 
k')ri]/m*. Expressions similar to Eqs.ljHEl '^''^^ be ob- 
tained for the advanced functions G^,^^{lu) = Gk/k(w)^. 
The sum over k" in the second Eq. |(SJ) can be directly 
calculated. First, we decompose G^„ into a spin in- 
dependent scalar part and a spin dependent part which 
is proportional to hk" • cr. Due to the time inversion 
symmetry hk" = — h_k" the sum over k" on the spin 
dependent part is zero for an isotropic scattering am- 
plitude. For anisotropic amplitude, however, this sum 
is not identically 0. Nevertheless, the sum on the spin 
dependent part can be ignored either way in the follow- 
ing calculations, because it is proportional to the small 
parameter kkf/Ep ^ 1. Further, it is easily seen that 
only Im[G|^„ ] is important in this k" sum because the 
real part gives rise to a term that simply adds to C/kk' in 
the first line of Eq. jSJ), thus effectively renormalizing the 
Born scattering amplitude. The imaginary part can not 
be absorbed in such a way because it has opposite signs 
for the advanced and retarded Green functions. Taking 
into account that lo ~ Ep and assuming that h^p <C Ep 
we thus get 

5]f/kk"G;;^?^(cc')[/k"k' = -z^^i.5(k,k')e'('''-''>' 

k" 

5(k, k') = -^ / d^"/(k", k)/(k', k") , (6) 

m J 

where </>" is the angle of the vector k", with |k"| = kp. 
At k = k' the integral in © is equal to the scattering 
cross-section. 

Within the semiclassical theory we follow the well 
known method [ij, |lj| to calculate the configurational 
average of the Green functions product that enters into 
the Kubo's hnear response equation. Because of scatter- 
ing on a target impurity this product contains more than 
a pair of Green functions. As our leading approximation 
we take into account the so called ladder perturbation 




FIG. 1: Diagram for the spin density. Scattering of electrons 
by a target impurity is shown by the solid circles. Dashed lines 
denote the ladder series of particle scattering by the random 
potential, p, k, k' are electron momenta. 



series describing particle and spin diffusion processes. 
When averaging the Green function product, within this 
approximation, only pairs of retarded and advanced func- 
tions carrying close enough momenta should be chosen to 
become elements of the ladder series. This considerably 
simplifies calculations. Some of the representative dia- 
grams are shown at Fig. 1, where v denotes the velocity 
operator 



k^ 



m 



9hk • cr 

dki 



(7) 



The diffusion ladder renormalizes only the vertex asso- 
ciated with the electric field, while such a ladder, as we 
just explained, does not appear at the vertex attributed 
to the induced spin density. It is because in the ballis- 
tic range around the impurity, the momentum transfer 
|p — k| ^ \I{vpt), and thus the diffusion is not impor- 
tant. Finally, the density of spins oriented in z direction 
can be written as 



'^.(r) = 



E 

k.k'.p 



j(p-k): 



duj dnp{uj) 



27r duj 
Tr[Gt,^{co)a,,Gl^,{u)T{co,k')], 



(8) 



where the trace runs through the spin variables and 
np{uj) is the Fermi distribution function. The functions 
Gk'k ^''^ given by Eq. Q. In (jSJl, only terms up to the 
second order in C/k'k should be taken into account. Hence 
the highest order corrections are those shown at Fig.l, 
b) and c). At low temperatures only uj in close vicinity 
around Ep contributes to the integral in (jHJ. Therefore, 
below we set u = Ep. 

It is important that the vertex T{uj, k) in Eq.© is the 
same that enters into the spin-Hall conductance. On the 
other hand, as was shown in many publications [l5j, for 
linear in k SOI and /ik ^ Ep, a contribution to T(uj,\i.) 
from the spin dependent part of the velocity Q cancels 
with its spin independent part, while such a cancellation 
does not take place in case of nonlinear SOI ,16]. As 



a result of such cancellation, in a linear case T(a;, k) is 
reduced to the simple expression 



result, which is as a sum of all diagrams in Fig. la)-c), 



T(c^,k) = — kE. 

TO* 



(9) 



Let us consider the spin density © in the presence of 
the Rashba spin-orbit field h^ = aky,hy = —ak^. In 
the zeroth order in f/kk' the Green functions in © are 
given by the first term in (PJ . In this approximation and 
with T given by Q one can easily see that crz(r) = 0. 
On the other hand, the inplane spin polarization directed 
perpendicular to E is finite. This polarization is due to 
the electric orientation effect [l7| . In the first order with 
respect to C/kk' the z spin polarization is represented by 
Fig. la. Expressing C/kk' via the scattering amplitude, 
from (|3I5|I and IJHI) we obtain 



<>{r) 



k,p 

(a,G^"(ii;^)/(p, k)e*(P-'')^ + h.c.)] , (10) 



where R = r — r^. At kR,pR ^ 1 the 2D angular inte- 
gration in (|10|l can be performed by expansions around 
the saddle-points (pR/pR) = ±1 and (kR/kR) = ±1, 
that result in the asymptotic expansion of ^^(r) at a 
large distance from the impurity. The remaining inte- 
grals over p and k are dominated by contributions from 
the poles of Green functions Q . These poles are located 
at k,p — zLkp ± Lj^ ± il^^, where Lso = H/m*a 3> kj, 
is the characteristic spin-orbit length and I is the mean 
free path. In the ballistic range R < I one may sub- 
stitute the imaginary part il^^ of the poles by iS with 
(5 — > 0. Depending on combination of ± signs of the 
poles, the scattering amplitude entering into IjlUI) will 
coincide either with the forward scattering amplitude 
/(O) = /(fciT-R, fci^R), or with the backscattering am- 
plitude /(tt) = f{kpR,—kpR), where R = R/R is the 
unit vector directed to the observation point. Finally, we 
obtain from itTUI) and igj 



4^Hr) = ^ 



R 

Re[/(7r)e 



f dnn 



vrafc.^Hs^'"^" 



2^fc.fllgij^2 



R 



T- - (11) 



where v^ = ctE/to,* is the drift velocity. The unit vector 
rifl = hj, ff^/\hp. f^\. For Rashba interaction it is n^ = 
Rv^ n\ = -i?^. 

In a similar way and with the use of Eq. @ one can cal- 
culate the second order contribution to the spin density, 
that is represented in Fig. lb) and c). Assuming the elec- 
tric field to be applied in the x-direction, we get the final 



^z{r) = - 



m v^at 




(12) 



where atot and at are the total and transport scattering 
cross sections, respectively, and 9 is the angle between the 
vector R and the x-axis. In order to check our method we 
applied it to the calculation of the charge dipole, whence 
az is substituted by 1 in Eq. ((SJ. Ignoring SOI we ob- 
tained the same result as in Ref . JJj . 

The explicit shape of the spin cloud is clearly seen from 
Eq. (|12|l . It consists of a dipole, oriented perpendicular 
to the electric field, and a tripole. Similar to the Lan- 
dauer charge dipole distribution [ll|, the spin density 
contains both slowly varying and fast Friedel oscillation 
components. Important distinctions, however, are found 
in the asymptotic behaviors. First, unlike the charge 
density, whose slow asymptotic term is represented by 
monotonous R~^ dependence, the spin density oscillates 
with a period determined by the spin-orbit precession 
length ttLso- Second, at smaller distances R < Lgo, the 
polarization decreases as R~^. It should be noted that 
this asymptotic form can not be obtained by the method 
based on the conventional leading order asymptotic ex- 
pansion of the wave function, as it has been done in |ll| 
for the Landauer dipole. It is because in 2D geometry the 
corresponding scattered amplitude decreases as 1/VR. 
Accordingly, the probability density, which can be either 
the charge or the spin density, will be proportional to 
l/R, not l/i?2. 

When talking about asymptotic expression H12|) . one 
should not forget that it is valid only within the ballistic 
range R ^ I. At larger distances the ballistic part of the 
spin density decays as exp{—R/l). On the other hand, 
outside the ballistic range the spin diffusion becomes im- 
portant. Spin diffuses during the D'yakonov-Perel' |l8| 
spin relaxation time, up to the distance ^ Lgo- Hence, 
the spin diffusion must be taken into account at R ^ I, 
providing that the spin-orbit coupling is not too strong, 
so that Lso ^ I- In order to calculate the spin density 
in the diffusive range, the ladder diagrams renormalizing 
the left hand vertex in Fig. 1 should be taken into ac- 
count. An evident result to be expected in this case is 
that the diffusion spin cloud with the size ^ I will ap- 
pear in addition to Ea. l|12|) . Due to the spin relaxation, 
however, the spin density will decay exponentially at 
R 3> Lso- This behavior is in sharp contrast to the power 
law decreasing of the charge density [lOj. In the latter 
case, the long-range R~^ charge-density tails of many im- 
purities result in the macroscopic electric field which can 
be related to the electric potential difference at the sam- 



pie boundaries. This was the main idea by Landauer - 
to associate impurities with resistors which give rise to 
an overall potential drop at a given current. Similarly, 
one would try to formulate the spin-Hall effect in terms 
of the spin-Hall resistance and spin dependent chemical 
potential x(r), defined as iVFx(r) — J2i'^z{i^ — i"i)- In 
such a way an influence of disorder on SHE can be con- 
sidered from the microscopic point of view starting from 
individual impurities. For example, similar approach has 
been employed in a semiclassical analysis of the side jump 
contribution to the anomalous Hall effect [l9j. Return- 
ing to the spin potential x(i"); one can notice that due 
to exponential decay in space of the spin cloud, the well 
converging sum over impurities will produce, on average, 
a vanishing spin-Hall chemical potential everywhere, ex- 
cept for the R ~ L^o range near the sample boundary. 
No such spin accumulation, on the other hand, has been 
found near hard wall flanks of a 2D diffusive strip of 
2DEG with Rashba SOI j3|. At the same time the finite 
accumulation was calculated in |a] for other boundary 
conditions. Probably, this means that the spin density 
outside the ballistic range around an elastic scatterer is 
finite, but the combined spin density produced by many 
impurities will depend on the boundary conditions for 
spin diffusion. 

In conclusion, for a 2DEG with Rashba spin-orbit in- 
teraction we calculated the nonequilibrium spin polariza- 
tion induced by the intrinsic spin-Hall effect in the ballis- 
tic range around a spin independent scatterer. The angu- 
lar spatial distribution of the spin density is represented 
by a tripole and a dipole oriented perpendicular to the 
electric field. As a function of the distance from the scat- 
terer, the polarization shows the power law decay with 
oscillations, some terms oscillating relatively slowly, with 
the period ttLso, while other terms varying fast, with a 
period of Friedel oscillations. Noteworthy, that although 
the z-polarized spin-Hall current is zero in case of Rashba 
SOI, we found out that the z-component of the spin den- 
sity is not zero in the ballistic range. This agrees with 
finite spin accumulation near flanks of a ballistic strip '8l| . 
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